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Abstract 

Quantum phases of naturally-occurring systems exhibit distinctive collective phenomena 
as manifestation of their many-body correlations, in contrast to our persistent techno- 
logical challenge to engineer at will such strong correlations artificially. Here we show 
theoretically that quantum correlations exhibited in the two-dimensional (2D) valence 
bond solid phase of a quantum antiferromagnct, modeled by Affleck, Kennedy, Lieb, 
and Tasaki (AKLT) as a precursor of spin liquids and topological orders, are sufficiently 
complex yet structured enough to simulate universal quantum computation when every 
single spin can be measured individually. This unveils that an intrinsic complexity of 
naturally-occurring 2D quantum systems — which has been a long-standing challenge for 
traditional computers — could be tamed as a computationally valuable resource, even if 
we are limited not to create newly entanglement during computation. Our constructive 
protocol leverages a novel way to herald the correlations suitable for deterministic quan- 
tum computation through a random sampling, and may be extensible to other ground 
states of various 2D valence bond phases beyond the AKLT state. 

Keywords: Quantum computation, 2D valence bond solid (VBS) phase 



1. Introduction 

Daily stories about a quantum computer have earned in our mind its image as super 
postmodern technology based on an artificial full control of quantum systems in the 
highest possible precision. Most current approaches to implement a quantum computer 
are based on a bottom-up idea in that we intend to build it by combining key elementary 
objects, as well summarized by the celebrated DiVincenzo criteria [1]. The idea implies 
that it should be equally, i.e., reversibly, doable to both create and annihilate at will 
many-body correlations (or entanglement if its nonlocal character has to be emphasized) 
among many qubits artificially. However, despite various promising candidates for a qubit 
and remarkable experimental progresses in their physical implementation, this unitary 
control of many-body entanglement in a scalable fashion is believed to remain the hardest 
challenge. 
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On the other hand, in a historical perspective, humanity has strived to find and tame 
resources present naturahy on the earth, rather than achieving our goal from scratch by 
our immediate ability. Speaking of energy resources for instance, instead of establishing 
self-contained energy cycle in the first place, we have learned to take advantage of more 
and more elaborated natural resources, e.g., from wood and water to oil and nuclear 
(although their sustainable use is still our long-term project). Here, based on the fact 
that nature realizes various quantum phases as manifestation of underlying many-body 
entanglement, we suggest taking a complementary, top-down vision in that we attempt 
to tame a resource of suitably structured entanglement, which could either exist in nature 
or be simulated relatively naturally within our technology. 

A key point is how to live with our limited ability. Once a specific natural resource 
of structured many-body entanglement is provided, we are supposed to utilize only op- 
erations which just consume entanglement without its new creation, such as local mea- 
surements and local turning off of an interaction. In this regard, our approach can be 
best compared to hewing out an arbitrary quantum evolving system from the "carving" 
resource, in enhancing Feynman's original intuition [2^ about quantum simulation in a 
computationally universal manner. A perquisite is that we could be content with limited 
technology to prepare only one specific resource of structured many-body correlation, 
that is less demanding than the full unitary control over arbitrary correlation. Moreover, 
if such a resource is available as a stable ground state of a naturally-occurring system, 
that gives us several possibilities to prepare, such as dissipative coolings or adiabatic 
evolutions, in a similar way that nature presumably does. Our target of the naturally- 
occurring two-dimensional (2D) system is the valence bond solid (VBS) phase of spin 
I's on the 2D hexagonal lattice, modeled by Affleck, Kennedy, Lieb, and Tasaki (AKLT) 
[31 13] , which is widely recognized as a cornerstone in condensed matter physics. Their 
VBS construction of the ground state in terms of the distributed spin singlets (or the 
valence bonds) has become one of most ubiquitous insights in quantum magnetism as 
well as in high-Tc superconductivity, and leads to modern trends of spin liquids and 
topological orders. 

It turns out here that the 2D VBS phase, represented by the AKLT ground state, 
provides an ideal entanglement structure of quantum many-body systems that can be 
suitably tamed through our limited capability to the goal of universal quantum compu- 
tation. Our top-down vision is materialized conveniently in taking advantage of a con- 
ventional framework of measurement-based quantum computation (MQC). Its methods 
have been developed to "steer" quantum information through given many-body correla- 
tions using only a set of local measurements and classical communication, under which 
entanglement is just consumed without new creation. Later we extend that in a wider 
program to tame naturally-occurring many-body correlations. 

In the context of MQC, the so-called 2D cluster state [S] is the first and canonical 
instance of such an entangled state that pertains to a universal quantum computational 
capability when every single qubit (spin |) is measured individually and the outcomes of 
the measurements are communicated classically [Hl[7]- Remarkably, it was already noticed 
in Ref. [8] that MQC on the 2D cluster state utilizes a structure of entanglement which 
is analogous to that of the aforementioned VBS state. Following such an observation, 
the tensor network states, as a class of efficiently classically parameterizable states in 
extending the VBS construction, has been used in Refs. [SI [TU] to construct resource 
states of MQC, where it was indicated that a certain set of the local matrices or tensors 
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Figure 1: A scheme of quantum computation through measuring the correlations of the 2D AKLT state, 
a representative state of the 2D VBS phase of spin |'s on the 2D hexagonal lattice. After a random 
sampling which assigns every spin per site to one of three axes x, y, and z, the typical configuration of 
the outcomes enables us to choose the backbone structure (described by a shaded pinkish-orange region) 
along which quantum computation is deterministically simulated in terms of a quantum circuit. Our 
protocol harnesses a pair (depicted as a dotted bond of the hexagonal lattice) of neighboring sites where 
one is measured in a standard basis and the other is done in a complementary basis, to accommodate 
the desired structure of space-time along the region of the backbone. An emergence of the time is 
simulated if both two bits of information out of measurements per site are communicated to the same 
direction (as depicted as the double arrows), on the other hand, an emergence of the space is simulated 
if two bits of information are communicated to the opposite directions (as depicted as a pair of the single 
arrows pointing apart). The figure, reminiscent of the Feynman diagram for scattering of two elementary 
particles, illustrates a microscopic view of the Figure [2] The two-qubit CNOT gate is implemented in 
the middle region between two quantum logical wires running from the right to the left. 

which describe the correlations (cf. Eq. ^ in our case) can result in a quantum unital 
map through the single-site measurement. Notably, however, most known examples 
considered so far, including additionally those e.g., in Refs. pTHTS] . are constructed to 
have such a convenient yet artificial property — as often referred as one of peculiar 
properties of the correlations of the 2D cluster state — that it is possible to decouple 
deterministically (by measurements of only neighboring sites) a ID-chain structure that 
encodes the direction of a simulated time as a quantum logical wire of the quantum 
circuit model. This peculiarity is said to be artifact of another less realistic feature of the 
2D cluster state in that it cannot be the exact ground state of any two-body Hamiltonian 
of spin i 's [T7| UH] . Thus one cannot expect such convenience for the correlations of a 
genuine 2D ground state of a naturally-occurring spin system. 

The main result of our paper is summarized in the following (informal) theorem and 
illustrated in the Figure [T] As elaborated in the text, we introduce a novel way to 
herald the correlations suitable for deterministic quantum computation through a ran- 
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dom sampling, to tame for the first time the genuine 2D naturally-occurring correlation. 
Otherwise it has natural tendency to split an incoming information into two outgoing 
information because of certain symmetric nature of the three directions at every site of 
the 2D hexagonal lattice. This seems to be the reason why MQC on the 2D AKLT state 
has been an open question in a long time, although the AKLT state by the ID spin-1 
chain was shown in Ref. [11] to be capable of simulating a single quantum wire of MQC. 

Theorem. Universal quantum computation can be simulated through consuming mono- 
tonically entanglement provided as the 2D AKLT state \Q) (defined as the VBS state by a 
spin I per site and described as a tensor network state of Eq. ^ ) of the size proportional 
to the target quantum circuit size. To this end, we leverage single-site measurements of 
every individual spin |, a bounded amount of classical communication of measurement 
outcomes per site, and efficient classical side- computation. 

We note that upon completion of the work, another approach to suggest computa- 
tional usefulness of the 2D AKLT state is presented independently in Ref. |TH] by Wei, 
AfHeck, and Raussendorf recently, through transforming the 2D AKLT state into a 2D 
cluster state using a clever mastery of the graph states (an extension of the cluster state 
to a general graph). In comparison, our proof constitutes a more direct protocol to use 
the correlations of the 2D AKLT state to simulate straightforwardly the quantum circuit 
model, with much less resort to the known machinery of the cluster state. Our construc- 
tion of the protocol also respects more explicitly the physics and topological nature of 
the 2D VBS phase, such as the edge states at the boundary and the (widely-believed) 
energy gap at the bulk. It clarifies their operational usage in quantum information in 
that the former is used to process the logical information of quantum computation in the 
degenerate ground subspace and the latter could contribute to its protection, in provid- 
ing some robustness against local noise, as the ground-code version of MQC discussed in 
the Section |6l 



2. 2D valence bond solid ground state 

The 2D VBS phase can be modeled by a nearest-neighboring two-body Hamiltonian 
of the antiferromagnetic Heisenberg-type isotropic interaction (i.e. J > 0) [4} [2T]. 
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where Sfe is the spin-| irreducible representation of su(2) at the site fc, and the summa- 
tion is taken over all the nearest neighboring pairs (fc, k') of spin |'s on the 2D hexagonal 
lattice. The particular weights to the biquadratic and bicubic terms are chosen conven- 
tionally to be the projector onto the subspace of the total spin 3 for every pair of (fc, k'). 
However the 2D VBS phase itself is supposed to persist around this AKLT point without 
a fine tuning of these weights, in the same way as the ID case. It is important to mention 
that our 2D VBS phase should be distinguished from the 2D valence bond crystal (VBC) 
phase, since VBC is usually used to refer to the phase that consists of the valence bonds 
in a broader sense, namely including not only the VBS phase but also the dimer phase 
etc. However, there are considerable differences between VBS and the dimer phase, for 
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example, regarding the global nature of entanglement and the origin of the ground-state 
degeneracy (with an open boundary condition). 

The 2D AKLT ground state is such a VBS wavefunction that the symmetrization of 
three (virtual) spin ^ 's to represent a physical spin | per site is made on a collection of 
the singlet pairs ^(|1^) ® |0^) — |0^) (g) |1^)) of spin 2 's, where each singlet is distributed 
along every bond of the 2D hexagonal lattice. The construction can be visualized like 
in the Figure 3.2 of Ref. [3] for instance. It is straightforward, for our convenience, to 
describe it as a tensor network state via the celebrated Schwinger boson method |20] , 



1^) = E 
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(2) 



where a^^k') the site k (or fc') runs over and the trace is taken by 

the contraction of the tensors according to their locations on the 2D hexagonal lattice. 
The boundary condition is assumed to be open in that we are simply given a finite bulk 
portion of the lattice, and, the boundary tensor B is set to be the identity according 
the [Appendix A[ The tensors at the site with the T-shaped or _L-shaped bonds are 
found to be given by 



[f1 



^T[r] = ^(io^)(ii®(oi+z®(i^i), 

Arhf] - (-|1^)(01 ^ (11 + ^ ® (01) , 
^T[-f1 = |l^)(0l®(0^|, 



Aj.[f] = -^(-r)(i^|®|i^) + z®r)), 

A^[-f] = -^(|l^)(0^|®|0^)+Z®|l^)), 
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respectively. Here the first ket and bra correspond to degrees of freedom by the left 
and the right, and the second ket or bra corresponds to that by the up or the down, 
respectively. The Pauh matrices are defined as Z = |0"^)(0^| - |1'^)(1''|, X = |0'=)(1'=| -f- 
|1^)(0^|, and Y = iXZ with the imaginary unit i — %/— T. An element of the Pauli group, 
including the identity, will be denoted as T later. In this tensor network description, the 
effective spin ^ labeled by |0^) , should better be understood as a manifestation of 
the fractionalized degree of freedom, the edge state, emergent at the boundary across a 
single bond. Later it will turn out that these tensors are interpreted as the logical action 
on these emergent edge states (in other words, degenerate ground states), when the spin 
is measured in the direction by its argument. 

The set of tensors in terms of the other bases has exactly the same structure (up 
to a possible overall phase) because of the rotational symmetry. In defining |0/1^) = 
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where a = |i^:~5i~f and our notation, for example for the basis, is meant to 
describe the tensors obtained in replacing |0/1^) and (0/l^| into |0/1^) and (0/l^| as 
well as Z into X in Eqs. ([3| and Q. 

The 2D AKLT state inherits various characteristics from the ID VBS state. For 
example, its correlations measured in terms of the two-point function was shown to 
decay exponentially in lattice distance with the correlation length f — l/ln(3/2) « 2.47 
PT] . Furthermore, in very contrast with the dimer phase, it also realizes a fractionalized 
degree of freedom on every boundary, called the edge state, as mentioned above. A 
numerical calculation [22] of the entanglement entropy confirms a qualitative nature of 
edge states. On the other hand, the spectral gap to the excited states in the 2D AKLT 
model is widely believed to persist in the thermodynamical limit, but has yet to be 
proved. 



3. Insight to the MQC protocol 

We intend to simulate the quantum circuit model through measuring the correlations 
at every site. We call the part of the 2D hexagonal lattice sites that corresponds to 
the quantum circuit (consisting of the quantum logical wires running almost horizontally 
and their entangling gates described vertically) a backbone, as seen in the Figure [2) The 
degree of the backbone site refers to how many neighbors it has along the backbone. The 
degree-3 backbone sites are used at every junction of the horizontal logical wire with a 
vertical entangling gate, so that they are required only occasionally. 

A key insight to construct our protocol is that since the reduced density operator of 
every spin | per site is totally mixed and isotropic, i.e., the normalized identity projector 
I;, we are able to extract 2 bits of classical information by measurements per site. Then 
it is sensible that in stead of obtaining them at once, a part of the information, indeed 
log2 3 bits in our case, is first extracted and we adapt the next stage according to it. 
It might be surprising that the first measurement induces a kind of randomization, but 
intuitively speaking, this part is crucial to separate the original quantum correlation that 
intrinsically involves genuine 2D fluctuations into the classical correlation (or, a statis- 
tical sampling) that can be still efficiently handled by a classical side-processor and the 
"more rigid" quantum correlation suitable for deterministic quantum computation. A 
global statistical nature of the AKLT correlations through the first stage guarantees, in 
an analogous way with the classical percolation phenomenon, that an embedding of the 
backbone (i.e. the target quantum circuit) can be found in the typical configuration of 
a heralded, randomized distribution of entanglement. At the second stage which imple- 
ments quantum computation, the protocol is invented in such a way that the standard- 
basis measurement and complementary-basis one, both of which are defined later, are 
paired (as depicted by the dotted bonds in the Figure [T]). This paired processing is cru- 
cial to make the logical information bound in the domain of the backbone, somehow in 
a reminiscent way how the deoxyribonucleic acid (DNA) holds genetic information in its 
famous double strands. 
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Figure 2: (Left) Illustrated is a typical distribution of matched bonds depicted as thicker (dark-cyan and 
light-purple) bonds. Here the bond is called matched when the pair of neighboring sites is assigned to the 
same axis through the polarizing measurement at the first stage of the protocol. The region to simulate 
a quantum circuit by the complementary bases, called the backbone, is going to be chosen by classical 
side-computation, so as to guarantee its global topology is identical to that of the target quantum circuit 
(cf. the Section [sj. To this end, the backbone must skirt some "off-limits" configurations of matched- 
bond clusters that induce additional undesirable logical loops if all of them are involved in the backbone. 
Accordingly, it is sufficient that every quantum logical v^ire is chosen from any horizontally percolating 
path that circumvents only some clusters of matched bonds, for instance, colored in light purple. (Right) 
A possible embedding of the backbone, described as a shaded pinkish-orange region, has been identified 
by analyzing classically the occurrence of matched bonds . The "branching-out" matched bonds that 
are involved in the backbone and described by a shaded yellow region are the extended region where the 
logical information of quantum computation is processed in the second stage of the protocol. However, 
by having no additional closed loop, they would not disturb the desired logical action. Based on the 
statistical property originated from the genuine 2D nature of correlations of the AKLT state, the almost 
sure success of such an identification of the suitable backbone is guaranteed in an analogous way with the 
emergence of the macroscopic percolating cluster in the bond percolation phenomenon. The microscopic 
view near the CNOT gate is highlighted in the Figure [T] 

Summary of the MQC protocol 

Now we outline our MQC protocol, which consists of two stages, (i) The first stage 
is to apply a measurement {M^ , M''^ , M^} which polarizes randomly toward one of the 
three orthogonal axes at every site. We define a degenerate projection M'^ {fi = x, y, z) 
as 

M.^y^(|f)(f | + (6) 

The set of {M^,M^,M^} constitutes the positive operator value measure (POVM) by 
satisfying J2fi=x y z -^^'^^ -^'-^'^ = 1, so that it is a valid local measurement with three 
alternative, random outcomes /i. We must record the outcome /i^ at every site k and 
collect the location of such matched bonds that the pair of the axes for the neighboring 
sites k, k' coincides, namely = A'fc' • Based on the occurrence of matched bonds 
(which need additional care in their use), we are able to determine the backbone by 
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efRcient classical side-computation in circumventing some rare "ofF-limits" configurations 
of matched bonds defined later. 

(ii) The second stage carries actual quantum computation, using further projective 
measurements at every site and feedforward of their outcomes. Once the backbone is 
identified, the computation is deterministic in a very similar way with MQC on the 2D 
cluster state. 



4. Widgets for the logical gates 

A basic idea of the widgets to construct the logical gates from the set of tensors is 
outlined first. Suppose the bond is not matched for a pair of the nearest neighboring sites, 
then either of the sites can be used as the part of the backbone, by treating the other as 
the non-backbone site. Let the assignment of the first-stage polarizing measurement at 
the backbone site be and the other be n), and, as an illustration, assume fx = z dX 
the T site and v = x dX the _L site. It turns out that, in measuring the non-backbone site 
in the basis of{(| |,(— | |} (since the state is already in this subspace by the preceding 
M^), we can always get a unitary action at the backbone site along the other remaining 
two bonds, regardless of which outcome has occurred. This is because conditioned on the 
specification of the tensor of Eqs. (4[5 ) to be, for instance, |0^) ® |0^) by the outcome 



(1^1 at the non-backbone site, the tensors at the backbone site are essentially fixed as 

iT[ii = r)(ii®(iio-) = -i3io^)(i^i, 

(7) 

iThi1 = |i^)(oi®(oio-) = ^|i^)ri. 

The overall numerical factor implies the fact that the outcome (|^| would occur with 

probability and of course the other outcome would do with the same probability i. 
Indeed, the measurements used for the unitary actions are unbiased in their output prob- 
ability. Here after, however, we do not describe the overall factor explicitly, since it is 
irrelevant to the statistics of the logical output of quantum computation. It can be readily 
seen that we can get a unitary gate by mixing these two tensors in such a way that for in- 
stance, measuring this backbone site in the basis {^((| | + (— I |)} 

provides X(= iril''] +AThf or XZ{^ -At[|''] +iT[-i'']) respectively Both uni- 
tary actions should be interpreted as the logical identity since the difference by a Pauli 
operator can be readily incorporated, as explained later, by the adaptation of the follow- 
ing measurement basis so that the Pauli operators are treated as the byproduct operator 
T. In more general, measuring in a basis on the plane spanned by the aforementioned 
two vectors is found to provide the logical rotation R'{9) = |0^)(0^| -I- e'^|l^)(l^| by an 
arbitrary angle 9 along the axis z. We call any of such bases complementary to the stan- 
dard basis {(f^l I ("I'^ll given the axis /i at the site. The choice of the complementary 
basis depends on the axis v of the associated non-backbone site. Accordingly, we may 
denote the complementary basis in an abstract way as {(7^'''|}, which will be explicitly 
defined later. 

Such a construction of a "mutually-unbiased" input to the backbone is the key tech- 
nique provided by the pair of the standard-basis and complementary-basis measurements. 
Indeed, this basic idea is shown to be extensible with an additional care even when the 
bond is matched. As it will get clearer, the set of sites measured in the complementary 



bases, which covers strictly the backbone, corresponds to the region where logical infor- 
mation of quantum computation is processed. Our construction enables us to decouple 
the logical output probability of quantum computation from the rest of the sites mea- 
sured in the standard-basis measurements (conditioned on the configuration of the axis 
per site after the first stage of the protocol). A suitable identification of the backbone is 
addressed later in the Section[5]as a global nature of computational capability, supported 
by a prescription to the matched bonds in the [Appendix B| 

Single- qubit logical gates 

Any single-qubit logical gate in SU (2) can be implemented by the sequence of the 
rotations R'^{9) — |0'^)(0''| -I- e'^|l^)(l'^| along two independent axes, using the Euler 
angles 6''s. Here we take jj. = z and x without loss of generality, and describe the detailed 
protocol for R^{9). 

The (9) can be only attempted if the axis /x of the backbone site is assigned to be 
z after the polarizing measurement at the first stage. Otherwise we teleport the logical 
information along the backbone by implementing the logical identity by measuring the 
site in the complementary basis with the fiducial angle (i.e., 9 = 0) until we reach the 
backbone site with /i = z. Let the axis of the associated non-backbone site be v, and the 
outcome when measured in the standard basis denoted as c € {0, 1} in corresponding to 
((— |. The complementary basis for the given axis /i = z is defined (at both the T 
and _L sites) as, \i v = x, 

- + (-1)'^-^') ((11 + (-11) + (1 - (-1)'^'') (- (f I + (-ri)] ' 

(8) 

where h S {0, 1} corresponds to each outcome of the orthogonal projective measurement 
with the angle 9 or + it respectively, and li v = y, 

= ^ [(1 + (-1)^'=^') H (f I + (-r I) + (1 - (-i)'e^') (f I + (-r I)] ■ 

(9) 

In all the cases, the logical action is given by 

^i[a](7"l''(6l)|M^|a) = XZ^®''R'{e), (10) 

where a partial contraction has been made, in the same way as in Eq. ([?]), with the rank-1 
tensor from the associated non-backbone site, using the notation 

A-j[a]=A^[a]\d'), A^[a\ = {c® V\ A^[a]. (11) 

This is the desired rotation up to a Pauli byproduct T = XZ^®'^ . Note that the choice 
of the various complementary bases only depends on the axis v and not on the outcome 
c at the non-backbone site, so that once the axis of every site has been assigned and 
the backbone structure has been determined, the backbone site and its associated non- 
backbone site can be measured in a parallel way. Likewise, the counterpart for R^{9) 
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Table 1: The indices (a^,a^) of the byproduct operator T = X°' for the complementary basis 

measurement {(7^'"!} given the axis fj,, where b is the outcome of the site measured in this basis, and c 
is the outcome of the associated site (connected to the former directly by an unmatched bond) measured 
in the standard basis of the axis u . 



can be obtained by 

(7^'^(^)i = [(1 + (-1)'^^') ((11 + + (1 - ("1)'^'') (- (i1 + (-11)] ' 

(12) 

= 273 [(1 + (-1)'^^'') (11 + (-11) + (1 - (-1)'^^') H (11 + (-11)] ■ 

(13) 

The measurement in these complementary bases result in TR^{9) with T ZX^®'^. 
Although it is straightforward to define (7^'''(0)| for R^{d), it is theoretically sufficient 
that the sites labelled by /x = y are only used for the logical identity (with 9 = 0), whose 
dependence to the outcomes is T = ^fcffic^&ffic®! r^j^^ dependence of the byproduct 
operator T on measurement outcomes is summarized in the Table [l] 

In the above, we have assumed that the (degree-2) backbone site has in tow its 
neighboring non-backbone site which is measured in the standard basis along v and 
provides the outcome c. If the backbone site is connected via the matched bonds to 
its associated non-backbone site, then the intermediate sites should be measured in the 



complementary basis too, and a correction to {c") (or (c ® l''!) of A[a] defined by Eq. ( 11 ) 
is required. Roughly speaking, the non-backbone sites connected by the "branching-out" 
matched bonds (colored in yellow in the Figure [2]) need to be treated in the same way 
as the backbone sites, so that we have to sum up all their byproduct operators. The 
motivation as well as technical detail of this prescription to the matched bonds is found 
in the [Appendix B| 

4-2. Two-qubit logical gate: CNOT 

Together with the arbitrary single-qubit logical gates on every logical wire, it is widely 
known that any entangling two-qubit logical gate is sufficient to achieve the universality 
of quantum computation when the quantum circuit is composed from the set of these 
gates. Here we take the target two-qubit logical gate to be the controlled-not (CNOT) 
gate, \0'){0'\^l + \V){V\(g)X. 

The two-qubit logical gate requires connecting a pair of the degree-3 backbone sites 
with the T and _L bonds for each. In particular, the CNOT uses the degree-3 T site 
with the axis = z and the degree-3 _L site with fi2 = x. One could in principle use 
a matched-bond cluster with the same axis, instead of the single degree-3 site. Observe 
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that 



At[7"I"(0)] = ^Tt (1 ® (0^1 + Z $5 (1^1) , 

(14) 

^±[7"''(0)] = j^T^ (1 ® |0^) + X ® |1^)) , 

where the byproducts are Tj- = XZ^^ (g) 1 and T_l = X''^®!^® Z, and the conditiomng 
axes v\ , vi of the complementary bases can be chosen freely, so that we have set v\ — x 
and V2 — z for convenience. So, when the vertical direction (the second degree of freedom 



in the tensor structure of Eq. (14)) is contracted through the sequence of the logical 
identities in the same way as in the part of a quantum logical wire, it is readily seen that 
this provides the desired CNOT (=\{\®\ + \® X + Z ®\ - Z ® X)). 

Here we need to calculate the total byproduct operator from this intermediate ver- 
tical part, namely T| = x'^i'"''^ Z^k^-k where the summation is taken over the degree-2 
backbone sites symbolically labeled by k. Needless to say, if there are, in between, the 
branches by the matched bonds, we have to take their contribution into account according 



to the prescription of the Appendix B in order to determine (af, af.) at every junction to 
the branch. Thus, T| essentially amounts to the whole accumulation of the byproducts 
from the vertical part. By absorbing Z^^ '^1®^ (which has included the additional one 
from T_l) into (0/l^| of At, as weh as into \0/r) of A±, we can reduce them 

to the byproduct operators along two logical wires. The whole logical action is TCNOT, 
where T = XZ'''^®^" <®'^ (g) X^^®^'' <®'^Z. 



4-.3. Initialization and readout 

The initialization and readout (in the sense of the quantum circuit model) of the 
logical wire can be simulated as well. Suppose, without loss of generality, these are 
always made in the |0/1^) basis. All we have to do is to find the degree-2 backbone site 
with /i = z on the rightmost (for the initialization) and on the leftmost (for the readout) 
of every logical wire, and to measure it in the standard basis {(| |}-As already 

described through the notation for the standard-basis measurement, according to Eqs. (|3j 
|4]), the outcome c must be interpreted as |c^) or (c ® V\. Note, in particular regarding the 
readout, that the outcome c at the leftmost boundary site is not immediately the logical 
outcome of the computation. The latter is evaluated only together with the byproduct 
operator at the end of the computation. 

4-. 4- Adaptation based on byproduct operators 

Since the byproduct operator T stays in the Pauli group, we can use the same machin- 
ery as the cluster-state MQC, to deal with the randomness of the measurement outcome. 
The key idea is to postpone the effect of T by adapting the following angles of the logical 
rotations. As an illustration, suppose we wish to apply the sequence of R^{9^) followed 
by R^{9^) up to some byproduct operator T, namely TR^{0^)R^{6^). Since 

(x^'®^'zi?^(r))(xz''®^i?^(r)) ^ x'''®'='®iz^®=®ii?^((-i)''®'=6i^)i?"(r), (15) 

we realize that if we adapt the second angle 6^ to be (—1)" 9^ , based on the byproduct 
index a" = & © c of the first measurement, we can always apply the desired sequence 
of the rotations regardless of the outcomes b and c. In general, we have to adapt the 
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angle of the next logical rotation, based on the current byproduct operator T updated 
according to the simulated time direction of the quantum circuit. Note that such an 
adaptation is required only if the rotation angle is not fiducial, so that the part by 
the logical identities does not need the adaption and can be implemented parallely in 
principle. Since this adaptation is now a widely-known machinery in MQC, we leave the 
details to the literatures, e.g., Ref. [7] (although in the case of the cluster state there is 
an additional automatic logical operation by the Hadamard matrix ^ {X + Z) at every 
single step of computation). 



5. Global nature of emergent computational capability 

5.1. Identification of the backbone 

In this section, we show that the backbone structure, which should be intact from 
rare "off-limits" configurations of matched bonds, can be identified efficiently in analyzing 
their occurrence by classical side-computation. The analysis not only guarantees that 
the widgets for the logical gates in the Section |4] can be composed consistently, but also 
suggests how the computational capability as a whole is linked to the emergence of a 
macroscopic feature of the many-body system. 

Recall that, at the first stage of the protocol, every outcome of the polarizing measure- 
ment {M^} = X, y, z) occurs randomly at every site, i.e., with the equal probability ^. 
Thus, in neglecting two-point correlation functions exponentially decaying in distance as 
a first approximation, we can imagine that the bond between a pair of the neighboring 
sites is unmatched with probability | and matched with probability |. Indeed, the effect 
of two-point correlations is very limited, and it seems to modify the probability of the 
unmatched bond by not more than 1 percent. Accordingly, we could visualize in our 
mind a typical configuration of unmatched bonds, as presented in the Figure [21 by the 
proximity to the bond percolation model with an occupation probability p — ^. As a 
reference, it is widely known in the bond percolation model that, when p is larger than 
the critical value Pc = ^ — 2sin(j|) « 0.652 . . . [53] for the 2D hexagonal lattice, there 
exists almost surely (i.e., with a probability close to the unity exponentially in the lattice 
size) a single macroscopic cluster of the occupied bonds, spanning the whole lattice. 



Importantly, it is shown (in detail in the Appendix B ) that not only unmatched 
bonds, which are the majority, but also most of matched bonds are indeed available as 
a part of the backbone under our prescription. The key idea is a "topological encoding" 
of the backbone structure in that in most cases the matched-bond cluster touching the 
backbone and measured in its complementary basis can be made to constitute a micro- 
circuit whose output to the backbone is mutually-unbiased as if the whole cluster were 
a single unmatched bond. In other words, the matched-bond cluster can be effectively 
renormalized to the single logical site directly on the backbone (which is called the "root" 



site in the Appendix B ), regardless of the size and shape of the cluster. This prescription 
clarifies that the backbone should be identified in order that the region measured in the 
complementary bases (namely, the backbone and all the matched-bond cluster touching 
it) does not induce any additional closed loop which is absent in the target quantum 
circuit. 

For a practical algorithmic implementation, we could leverage a variety of polynomial- 
time classical algorithms (see e.g., Ref. [23]) developed through the study of percolation 
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models. For instance, using an existing method to detect the locations of matched-bond 
clusters locally from the records of the polarizing measurement outcome at every site, we 
can efficiently mark potential off-limits configurations that two matched-bond clusters 
with different axes are nearest neighboring via more than one unmatched bond (see 



the Appendix B for their precise notion related to the existence of closed loops). This 
is because measuring both clusters in complementary bases must result in an undesired 
closed loop inevitably. So, one could consider one of such clusters to be effectively 
"unoccupied" bonds (which are marked by thicker, light-purple bonds in the Figure [2]), 
and choose the backbone from the remaining (unmatched and matched) bonds. 

Thus, our situation is much more favorable to find every horizontally percolating path 
corresponding to a quantum logical wire, compared to the case where such a path should 
be found from the cluster of unmatched bonds only. The ratio of "occupied" bonds is 
much larger than | effectively, and our macroscopic cluster of the occupied bonds can 
indeed offer abundant choices. For this purpose, it is straightforward to apply to our con- 
text the same efficient classical algorithm to find multiple percolating paths, despite that 
the unoccupied bonds are not generated according to the standard percolation model. 
That may involve, in practice, some customization such as a preference to the sites with 
certain axes (depending on the choice of the elementary gates) and optimization about 
physical overheads, but we do not intend to pursue concrete algorithms further. Note 
that the use of a percolation phenomenon has been previously considered in the context 
of MQC but in rather different scenarios [25l [26] where the preparation of a cluster state 
suffers some locally probabilistic errors. For example, Ref. [2(3 has presented a detailed 
algorithm to extract explicitly multiple percolating paths for a 2D cluster state. Such 
machinery to the cluster state is also used in another approach to indicate computational 
usefulness of the 2D AKLT state in Ref. [T^. 

Only a practically notable constraint by the matched bonds is on the spatial overhead 
of the backbone, in particular, the spacing between every neighboring pair of quantum 
logical wires running horizontally. Two quantum logical wires are not supposed to enter a 
single matched-bond cluster, since the latter behaves as the single logical site effectively. 
According to the approximation by the percolation model, it is expected that the size of 
the largest matched-bond cluster scales logarithmically in the lattice size but it occurs 
with an exponentially small probability. Since we would likely resort to quantum error 
correction to the goal of a scalable fault-torelant quantum computer in the end, we can 
set a cutoff within its constant error threshold, so as to disregard the rare large clusters 
and to deal with them in the similar way with other errors. It means that a constant 
spacing between two neighboring logical wires would be sufficient. 

5.2. Emergence of a simulated space-time 

It is remarkable to see how the distinction between (simulated) space and time 
emerges in our model. Our backbone structure after the first stage does not distin- 
guish them yet, but we can realize, by a careful examination of the protocol, it originates 
from the way the byproduct indices (a^, a^), namely two bits of classical information pro- 
cessed at each site, are communicated to the neighbors. Time emerges when both two 
indices are forwarded in the same (and roughly horizontal here) direction, while space 
emerges when two are sent in the opposite (and roughly vertical here) directions as is the 
case of CNOT, where or is sent downward or upward, respectively. Note also that 
since the spatial-entangling gate, like CNOT, is simulated only using the measurements 
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with the fiducial angle and thus need not be adapted, there is indeed no "time-ordering" 
along the simulated spatial direction. 

6. Discussion 

6.1. Ground-code MQC 

A feature that our resource is available as the ground state can be leveraged during 
the computational process in addition to the stage of its preparation. As originally 
proposed in Ref. [TT] as the ground-code version of MQC using coupled ID AKLT states, 
if we are capable of adiabatically turning off a two-body interacting term only to the 
spin to be measured selectively — importantly that still does not have to create new 
entanglement — we could indeed make the Hamiltonian in the remaining bulk coexist 
without interfering the MQC protocol for quantum computation. Then, the Hamiltonian 
with a (conjectured) gap would be able to contribute to a passive protection of the logical 
information stored in the degenerate ground state, by providing some robustness against 
local noises, in compared to a scenario where the bulk Hamiltonian is absent. 

Here we need additional care to apply the idea of the ground code to the 2D AKLT 
state, since it was assumed that the spin |'s are first assigned to one of three axes through 
the polarizing measurement before the quantum computation starts. It turns out that 
it is still possible for the the classical algorithms described in the Section [5] to function 
if the sampling of the matched bonds are made over the block of a constant number of 
spins near the rightmost boundary of the unmeasured bulk, without identifying the full 
configuration of matched bonds deep inside the bulk. 

6. 2. Persistence of computational capability in the VBS phase 

Another significant merit to have the bulk Hamiltonian during the computation is 
that any ground state which belongs to the 2D VBS phase (with the rotational sym- 
metry as is present in the AKLT Hamiltonian) is conceived to be ubiquitously useful to 
simulate the same quantum computation through the aforementioned adiabatic turning 
off of the interaction at the boundary. Its detail will be analyzed elsewhere, but such 
persistence of computational capability over an entire quantum phase was unveiled in 
Ref. [57], in examining the rotationally-invariant ID VBS phase to which the ID AKLT 
state belongs. A key of the mechanism is indeed the edge states — fractionalized degrees 
of freedom emergent at the boundary of the unmeasured bulk part — that are commonly 
present in the VBS phase and carry the logical information of quantum computation in 
our context. At any point within the VBS phase, the quantum correlation between the 
spin to which the two-body interaction has been turned off and the emergent edge state 
at a new boundary of the bulk would be modified to exactly that of the AKLT state, 
since the AKLT Hamiltonian satisfies additionally a frustration-free property (that the 
global ground state also minimizes every summand of the Hamiltonian) among a param- 
eterized class of generally-frustrated Hamiltonians within the phase. That is how the 
primitive of the ground-code MQC, the adiabatic turning off of the interaction followed 
by the measurement of the freed spin, can ubiquitously unlock computational capability 
persistent over the VBS phase. 

In contrast, suppose our ability of the selective control of a single spin at the boundary 
is further limited only to its measurement and the bulk Hamiltonian is neither engineered 
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nor present during the computational process. We may still consider implementing phys- 
ically (namely probabilistically in compensation with a longer distance scale) a quantum 
computational renormalization, through which the correlation within the VBS phase is 
modified into that of a fixed point. This is envisioned based on Ref. [28 where a protocol 
within the ID VBS phase was constructed so as to set the AKLT point as its fixed point. 
An attempt to construct a 2D counterpart of the protocol may face the intrinsic difficulty 
for the traditional computer to analyze correlations in the 2D system, but recent various 
developments of the classical methods (cf. in Refs. PUH32) ) regarding the renormalization 
of 2D entanglement might be applicable to our context. 

6.3. Outlook towards physical implementations 

The persistence of such computational capability over the 2D VBS phase would be 
encouraging towards physical implementations of our scheme. Then a fine tuning to 
engineer a set of specific interaction strengths in the 2D AKLT Hamiltonian of Eq. ([!]) 
may not be necessary as far as the system can be set within the "quantum computational" 
phase in maintaining some key symmetries such as the rotational symmetry. 

The biggest difference from the ID case is the fact that the ordinary Heisenberg 
point (i.e. the ground state of the Hamiltonian without biquadratic and bicubic terms 
in Eq. ([T])) has been known to be Neel ordered ,4J and does not belong to the 2D VBS 
phase unfortunately, although the ID Heisenberg point is located within the ID VBS 
phase. That implies that we may need a more elaborated interaction than the simplest 
Heisenberg interaction or a use of geometrically frustrated quantum magnets, in order to 
set a ground state in the 2D VBS phase. In fact, the ID VBS phase has been identified 
experimentally in several natural chemical compounds such as CsNiCls, Y2BaNi05, and 
so-called NENP. The 2D VBS phase seems to be less studied both theoretically and 
experimentally, compared to other kinds of 2D VBC phases for now. A recent discovery of 
the absence of spontaneous magnetic ordering in the manganese oxide Bi3Mn40i2(N03) 
[55] , which materializes an antiferromagnet of spin | 's in terms of Mn^+ ions on the 2D 
hexagonal lattice, might pave the way towards this direction (e.g., Ref. [51]). 

Another possible realization of the 2D VBS phase is based on analog engineering of 
the Hamiltonians of spin lattice systems. A recent experimental development of quantum 
simulators in atomic, molecular, and optical physics is promising in that Mott insulating 
phases of ultracold Fermi spinor gases with a hyperfine manifold F — e.g., in terms of 
^Li or ^"^^Cs atoms trapped in the optical lattice, might be used to simulate the 2D AKLT 
Hamiltonian of spin | 's along an extension of the approach in Ref. [35; for example. 
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Appendix A. Boundary tensor in the 2D VBS state 

The boundary tensor B is defined as _B = J^^ T^, where i represents a nonlocal degree 
of freedom associated with every pair of the boundary sites which are supposed to be 
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closed under the periodic boundary condition, and thus 2\£\ represents the total number 
of the boundary sites. This term formally provides the 41^' -fold initial degeneracy to the 
ground state, but it turns out that any of such ground states (and actually even the 
mixed state by them through an argument similar to Ref. |27j ) is useful to our goal, 
because of the decoupling property to the backbone. Thus, here we set B to be the 
identity for simplicity. 

Appendix B. Renormalizing prescription to matched bonds 

We elaborate a technical prescription for matched bonds, and prove the following 
informal lemma. 

Lemma. The backbone is freely chosen as far as the global topology of the region measured 
in the complementary bases (namely, the backbone and all the matched-bond clusters 
touching to it) is equivalent to the target quantum circuit, regarding the configuration of 
closed loops. 

The lemma can be also seen as our resolution to the genuine 2D fluctuations present 
in naturally-occurring systems. They are dealt, by adjusting the encoding of each single 
logical site to be variable according to the size and shape of every matched-bond cluster. 

Appendix B.l. ID-like chain 

Let us consider first a ID-like chain of the matched bonds, branching out from a 
degree-2 site on the backbone of the axis /i. Here we call such a site as the root site of 
the matched-bond cluster. The sequence does not bifurcate further, and thus is expected 
to terminate in the end at a non-backbone site after a chain of matched bonds. Let 
us for now label these sites along the matched-bond chain by k. As said, all the non- 
backbone sites along the chain of the matched bonds are measured in the complementary 
bases {(7''''^'' |}, which depend on the axis of each associated site connected by the 
unmatched bond in the third direction orthogonal to the matched-bond chain. At the 
end of the matched-bond chain, the terminating site must have two unmatched bonds, 
and its two associated sites are measured in the standard basis of each. That provides a 
pair \c'^) and jc'" ) to the tensors at the terminating site, so that either one, here jc'" ), 
can be used to assign A[a\ to the terminating site as usual, and the other \c^) plays a 
role of the "initial" condition to the matched-bond chain. In the exactly same way as the 
chain of the degree-2 backbone sites, this matched-bond chain behaves as the sequence 
of the logical identities acting on the initial state jc"). At the root site to the backbone, 
the total logical effect is equivalent to input \c'^) given by 



matched 



(B.l) 



k 



where u = v and 




(B.2) 
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The indices (0^,0^) of the byproduct operator at every site k, measured in the comple- 
mentary basis along the matched-bond chain, are determined according to the rule at 
the Table [T] (in incorporating the contributions from associated sites measured in the 
standard basis), (c© l^j is determined exactly in the same formula. The key observa- 
tion is that the matched-bond chain all measured in the complementary basis provides 
a "mutually-unbiased" input state in an axis v different from the axis fi of the matched 
bonds, and thus can play a comparable role to a single unmatched bond. 

Appendix B.2. Locally tree structure 

In general, the branching-out cluster of the matched bonds may have bifurcations like 
a tree locally. Flowing backward from every leaf of the tree (or, the terminating site of 
a matched-bond chain) according to the aforementioned procedure, we can prepare the 
mutually-unbiased input at each bifurcation site. Thus, we concatenate this procedure 
several times until it reaches the root site of the matched-bond cluster located on the 
backbone. 



Appendix B.3. Closed loop 

For completeness, we remark that closed loops may exist within the matched-bond 
cluster, branching out from the backbone. It essentially corresponds to the counterpart of 



Eq. (B.l I in a closed boundary condition. Suppose a closed loop, for example a hexagon. 



exists when we are applying the aforementioned procedure from every leaf of the cluster. 
Each site in the closed loop, except the the one which is closest to the root site of the 
cluster and called as the zeroth here, is now assumed to be associated effectively with the 
mutually-unbiased input using either an unmatched bond or our prescription for matched 
bonds. Then, the contribution by the closed loop is 



tr loop 



matched 



(B.3) 



where the axis uq can be chosen to be either axis different from /i, and note that Aq is 
the original, three-directional tensor. For example, if /i = z and the zeroth is a T site, 



Eq. (B.3 1 is equivalent to 



c=0,l 



(B.4) 



by denoting the zeroth-site bare outcome as 69. Remarkably, careful examination of the 
Table [T] clarifies that the total byproduct operator along the loop could accumula te up 
to the identity operator (with a regular weight ^), and the contribution of Eq. ( |B.3[ ) 
is always non- vanishing. In this example, 1 © = because the number of sites 

constituting a loop must be even and every a% is always 1 according to the Table 
Thus, the closed loop of matched bonds provides likewise a mutually- unbiased input |c' 
or (c© 1"! towards the root site, where v = v^) and 



60 © Efc al 
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if /i = 
ii H = z, 



(B.5) 



in analogy to Eq. (B.2) 



This analysis suggests that the closed loop itself can be directly a part of the back- 
bone, namely can share one or more bonds with the backbone. Phrased differently, the 
matched-bond loop (of the same axis) does not disturb the global topology of the back- 
bone. This can be confirmed by inserting, instead of Aq, two original three-directional 
tensors, e.g. Ai^ and Aout at a T and a _L site respectively with fi = z and vq = x or y, 
in Eq. (B.3). Note that the two are not necessarily the nearest neighbors, nor necessarily 



T and _L sites in general. 



tr 



loop 



matched 

Aout [7"'"° (0)] n ^[7''"'=(0)]A„[7"''^°(0)] 



c,c'=0,l 



\c"'°)(ci 



where in the last equality we have used the non- vanishing condition of the trace, c© c' = 
b-f ® b± (B J2k'^k (i-^-' parity of c and c' is fixed), based on the property that the 
exponent of the Pauli X operator is identically zero once again. Thus, the contribution 
by the loop is a power of the Pauli fi operator depending on the measurement outcomes 
in an extended manner of the Table [T] It is also straightforward to see that we could get 
further R^{6) by this matched-bond cluster, by replacing one of fiducial rotations by 6, 
as usual. 

In a nutshell, as mentioned in the main text, the cluster of matched bonds can be 
effectively treated as a single logical site under our renormalizing prescription, regardless 
of its size and shape. 



Appendix B.4- Off-limits configuration 

On the contrary, suppose a closed loop of sites measured in the complementary bases 
is introduced by a situation that the backbone touches sequentially two kinds of matched- 
bond clusters (with different axes, say ^ = x and z) nearest neighboring via two or more 
unmatched bonds. Then the contribution of the induced loop does change the global 
topology of the backbone. As an illustration, suppose the induced loop consists of the 
matched-bond cluster of the z axis and a single site of the x axis, all of which are 



measured in their own complementary bases. By replacing in Eq. (B.6) the axis z at A 



by X (implicitly there are two unmatched bonds in the loop), we see that 



f r loop 



matched 



Aout [7"'"" (0)] n Afc[7''''H0)]Ai„[7"l''"(0)] 



tr |^X''^®'=®^'=''^Z''^®'='®^fc°fc 

c,c'=0,l 



(B.7) 



The resulting map is not unitary anymore in contrast to Eq. (B.6). Such an off-limits 



configuration would be readily circumvented through local analysis of the configuration 
of matched-bond clusters. In the Figure [2] we have illustrated marking (by light purple) 
the off-limits configurations, where each cluster is paired with another cluster in such a 
manner that both should not be measured in the complementary bases. 
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